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ABSTRACT 


The  oblique  motion  of  a  circular  cylinder  through  an  inviscid  and  incompressible  fluid, 
conveyed  by  a  uniform  flow  at  infinity,  in  the  vicinity  of  another  cylinder  fixed  in  space  is 
considered.  In  a  relative  polar  coordinate  system  moving  with  the  stream,  the  kinetic  energy  of  the 
fluid  is  expressed  as  a  function  of  six  added  masses  due  to  motions  parallel  and  perpendicular  to 
the  line  joining  the  centers  of  the  cylinder  pair.  The  Lagrange's  equations  of  motion  are  then 
integrated  for  the  trajectories  of  the  moving  cylinder.  In  order  to  evaluate  the  added  masses  and 
their  derivatives  with  respect  to  the  separation  distance  between  the  cylinders  in  terms  of  the 
hydrodynamic  singularities,  the  method  of  successive  images,  initiated  by  HicksHl,  and  the 
Taylor’s  added-mass  formula  are  applied,  and  analytic  solutions  in  closed  form  are  obtained 
thereafter.  The  dynamic  behavior  of  a  drifting  body  in  the  close  proximity  of  a  fixed  one  is 
investigated  by  considering  the  limiting  values  of  the  fluid  kinetic  energy  and  the  interaction  forces 
on  each  body.  The  reliability  of  the  numerical  approximation  of  added  masses  and  their  derivatives 
is  also  discussed  in  the  present  study.  The  integral  equations,  in  terms  of  surface  source 
distributions  and  their  derivatives  on  both  circles,  are  carefully  modified  for  obtaining  accurate 
numerical  solutions. 
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I.  INTRODUCTION 


The  problem  presented  by  motions  of  two  or  multiple  solid  bodies  in  an  inviscid  and 
incompressible  fluid  anoears  in  many  research  fields.  For  example,  the  determination  of 
hydrodynamic  interactions  between  an  offshore  structure  and  floating  ice  masses  is  one  of  the 
practical  problems  in  the  oil  industry.  The  prediction  of  the  amount  of  liquid  droplets,  which  are 
conveyed  by  the  air  flow  and  collected  by  a  solid  surface,  is  important  in  the  atmospheric  icing 
research.  In  all  cases  concerned,  the  Reynolds  number,  which  is  the  ratio  of  the  inertia  effect  to  the 
viscous  effect,  is  fairly  large  such  that  the  potential-flow  theory  is  reasonably  applicable. 

One  and  half  centuries  ago,  Hicks^l  first  considered  motions  of  two  circular  cylinders 
moving  in  any  manner  in  an  unbounded  fluid.  He  studied  the  velocity  potential  due  to  the 
distribution  of  unit  sources  spread  over  each  surface  which  was  first  assumed  to  be  stationary,  and 
then  found  the  velocity  potential  due  to  the  motion  of  two  cylinders  by  proportionating  the 
magnitude  of  each  source  to  the  normal  motion  of  the  surfaces.  His  method  may  be  too 
complicated  to  be  applied  in  such  cases  where  the  main  concern  is  the  motion  of  solids  only.  A 
general  framework  for  the  motion  of  solids  in  an  inviscid  fluid  was  presented  by  Lamb^f  He 
applied  the  Lagrange's  equations  of  motion  in  generalized  coordinates  to  the  hydrodynamic 
interaction  problem  and  related  the  fluid  inertia  to  the  equations  of  motion  by  virtue  of  the  kinetic 
energy  of  the  fluid.  As  an  example,  he  considered  one  sphere  approaching  another  along  their 
centerline  in  a  uniform  flow.  Yamamoto^]  and  Isaacson  &  Cheungt4)  obtained  analytical  and 
numerical  solutions  for  forces  on  a  pair  of  cylinders.  Yamamoto^]  represented  the  complex 
potential  field  in  terms  of  an  infinite  series  of  doublets  by  the  method  of  successive  images,  which 
was  first  introduced  by  Hickst5!  for  the  motion  of  a  pair  of  spheres.  Yamamoto  also  derived  a 
formula  for  forces,  with  some  mistakes^,  on  each  cylinder  based  on  the  Blasius  theorem. 
Isaacson  &  Cheungf4l  did  the  same  work  numerically  and  their  results  agreed  very  well. 
However,  their  derivations  about  forces  on  solids  may  hardly  be  generalized  to  the  three- 
dimensional  problems  because  of  the  invalidity  of  the  Blasius  theorem.  A  better  way  of  describing 
the  motion  of  solids  would  be  to  use  the  Lagrange’s  equations  of  motion  and  to  evaluate  the 
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hydrodynamic  interactions,  in  terms  of  added  masses,  by  means  of  the  generalized  Taylor’s  added- 
mass  formula.  The  application  of  the  method  of  successive  images  is  limited  to  certain  symmetric 
geometries  such  as  a  pair  of  spheres  or  circular  cylinders,  due  to  its  dependence  on  internal  isolated 
singularities  which  cannot  be  found  in  general.  Landweber  and  Chwangl7^  developed  a  general 
formula  and  a  numerical  model,  based  on  the  boundary  integral  method,  to  obtain  the  added 
masses  and  added  moments  of  inertia  by  solving  the  unknown  source  distributions  on  solid 
surfaces.  When  two  bodies  are  not  quite  close  to  each  other,  their  method  gives  good  numerical 
results  on  added  masses  for  various  geometries.  However,  this  technique  depends  on  the 
simultaneous  solution  of  a  set  of  Fredholm  integral  equations  of  the  second  kind  which  may  cause 
certain  inaccuracy  when  the  bodies  are  very  close.  There  is  room  for  further  research  to  modify  the 
integral  equations  in  order  to  obtain  a  better  result  for  added  masses  and  their  derivatives  with 
respect  to  the  separation  distance  in  a  region  where  the  inviscid  and  irrotational  theory  itself  is 
valid. 

Motivated  by  the  lack  of  studies  on  the  oblique  motion  of  solids,  we  have  investigated  "he 
hydrodynamic  interactions  due  to  the  off-center  motion  of  a  circular-cylinder  pair  and  obtained 
analytical  and  numerical  solutions.  The  motion  of  each  body  is  assumed  to  be  purely  translational 
because  at  any  time  instant,  the  moment  due  to  hydrodynamic  interactions  is  zero.  The  kinetic 
energy  of  the  fluid  is  first  expressed  in  a  polar  coordinate  system  moving  with  the  stream  in  terms 
of  six  added  masses  of  the  bodies  due  to  their  motions  parallel  and  perpendicular  to  the  line  joining 
the  centers.  These  added  masses  are  in  turn  represented  in  a  relative  rectangular  coordinate  system 
as  a  function  of  the  separation  distance  and  of  the  angle  between  the  uniform  flow  and  the  polar 
axis.  By  relating  the  absolute  and  relative  frames  of  reference,  the  equations  of  motion,  in 
Lagrangian  form,  are  then  integrated  numerically  for  trajectories  of  the  drifting  body  in  space. 

In  Section  2,  we  are  concerned  with  the  Lagrange's  equations  of  motion  expressed  in  the 
relative  and  the  absolute  rectangular  coordinates.  The  transformation  of  added  masses  from 
rectangular  coordinates  to  polar  ones  is  also  considered  in  this  section.  Section  3  is  mainly  devoted 
to  the  exact,  closed- form  solutions  of  added  masses  and  their  derivatives  for  a  pair  of  cylinuas  due 
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to  centroidal  and  transversal  motions  respectively.  In  Section  4,  we  have  studied  the 
corresponding  numerical  solution  based  on  the  boundary  integral  method.  The  integral  equations 
for  the  source  distribution  on  each  surface  are  carefully  modified  and  the  solutions  are  compared 
with  the  exact  results.  The  trajectories  of  the  moving  body  under  various  conditions  are  given  and 
discussed  in  Section  5.  Finally,  conclusions  are  presented  in  Section  6. 

II.  EQUATIONS  OF  MOTION 

Consider  the  motion  of  two  circular  cylinders  in  an  unbounded  fluid  which  moves  with  a 
uniform  velocity  U0  in  the  x  direction  at  infinity.  In  an  absolute  frame  of  reference  fixed  in  space, 
body  1  with  center  oi  moves  with  velocities  u  and  v  in  the  x  and  y  directions  respectively,  and 
body  2  with  center  02  is  fixed  in  space.  Relative  to  a  moving  frame  of  reference,  in  which  the  fluid 
is  at  rest  at  infinity,  the  center  Oi,  located  at  (X[,yi),  moves  with  velocities  Uj  and  U3  along  the  x 
and  y  directions  respectively,  while  the  center  02,  located  at  (x2,o),  moves  with  velocity  U2  along 
the  x  direction  (see  Fig. la), 

U 1  =  u  -  U0,  U3  =  v,  U2  =  -U0.  (1) 

The  separation  distance  s  between  the  centers  is  given  by 

s2  =  c2  +  h2,  (2) 


where  c  =  .>q  -  x3  and  h  =  yi  . 

The  fluid  is  assumed  to  be  inviscid  and  incompressible  and  the  flow  irrotational.  Hence, 
there  exists  a  velocity  potential  0  which  may  be  expressed  as 

0  =  UiOi  +  U202  +  U303,  (3) 

where  Oi  is  the  unit  velocity  potential  corresponding  to  Uj  =  l,  U2  =  0  and  U3  =  0,  etc.  The 
kinetic  energy  T  of  the  fluid  can  then  be  expressed  as 
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(4) 


2T  =  A11uf+A22U2+A33U3+2A12U1U2+2A13U1U3+2A23U2U3, 
in  which  Ay's  (i,j  =  1,2,3)  are  added  masses  given  by 


3d, 

Ay  =  -  p  f  Oir^dS 
S’j  5n' 


Ay  -  Aji 


(no  sum  on  j), 


(5) 

(6) 


p  is  the  mass  density  of  the  fluid,  n3  =  ni,  S3  =  Si,  nj  denotes  distance  along  the  outward  normal 
to  the  jth  body,  and  the  integral  extends  over  its  surface  Sj. 

The  motion  of  body  1,  due  to  the  hydrodynamic  interaction,  can  be  determined  from  the 
Lagrange's  equations  of  motion. 


d  3T 
dt  aui 


+ 


(7a) 


M 


dU3  d  3T 

1  *  “'dtau3 


(7b) 


where  Mi  is  the  mass  of  body  1  and  t  is  the  time.  The  added  mass  Ay  depends  on  the  separation 
distances  c  and  h.  From  (2),  we  have 


9T  JT  3T=3T  dc  r  dh  _ 
3x,  dc  ’  ayi  ah  ’  dt  r  2’  dt 

From  equations  (4), (7),  and  (8),  we  have 


(8) 


Mi 


dUi 

dt 


l3Allu2  +  (aAu+13Aj?)u2_{aA13_iaA23)U2 
^  3c  1  3c  2  3c  2  3h  ^  3c  3 


dAn 

3c 


UiU2 


dAj_l 

3h 


U,U3+( 


3c 


3A 


12  ,  3A23mt  it  -  dUl  ■  *  dU?  ■  *  dU3x 
+  )U2U3  -  (An-^-  +  A  12-^r"  +  A13  -^7-) 


3h  3c 


dt 


dt 


(9a) 
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dU3  1  3AU  9A13  2  9A23  l  3A22  .  .  .2  1  9A33  2  9A33 

Ml— 7—  =  (o  -  -  -  )U  +( -  +7T -  )U  -7T - U  - UiU3 

dt  2  3h  3c  l  9c  2  ah  2  2  9h  3  9c 


3A;3  dA-n,  3Ap  ....  9A-*3  dU i  dUo  dU3 

+(  “  +  ^7^  )  k  ]U2  +—7^-  U2U3  -  (Ai3— +  A23~^p  +  A33  . 


dc  dc  3h 


dc 


From  (1),  we  have 


dU; 

dt 


=  0, 


dUi  du  dU3  dv 


dt  dt  ’  dt  dt 


=  7,7,  U3=v  . 


(9b) 


(10) 


By  means  of  (1)  and  (10).  we  may  simplify  equations  (9a, b)  as 


du  dv  1  9A|i  2  ^An  1  9A33  n  2  3Aii 


(M1+a„)^  +  a13^  =  -5^-( 


X 


dt  —  dt-2  3c  u  ah  2  ac  '  ah 


uv 


l3Au  1  3A22  3A12  2  ,3Au  9a12  9Ai3  9A23 

+()  — +  i—  +  — )U  +(-—  +  -— --^—--——)U0v, 
~  3c  1  dc  3c  o  3h  9h  9c  9c 


,u  .  A  xdv,  A  du  ^  1  3A ! !  9AI3  2  1  9A33  2  9A33 
(Mi  +  A33)  ^-  +  An  ^-  =  (  j  -7—  -  — — )u  -  o— — v  -  - — uv 


(11a) 


9h  9c 


9h  9c 


+  ln 


dA  1  j  1  3A22  9Ai2  2  ,9a  1 1  9 Aj~  3Ai3  9A23 


=)U  -( 


)UQu 


,  _  -)  ■»"  /  "  V  ~T  -  - 

“  dh  “  9h  9h  O  dh  9h  3c  9c 


(lib) 


The  added  masses  in  equations  (9a, b)  and  (1  la,b)  are  functions  of  c  and  h  and  evaluated  in  the 
rectangular  coordinates. 

For  a  pair  of  circular  cylinders,  we  may  represent  the  motion  of  body  1  in  a  polar 
coordinate  system  centered  at  (x2,  0), 

x  1  =  x2  +  s  cosy, 

y  1  =  s  siny,  (12) 
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and  decompose  each  added  mass  into  a  part  depending  on  the  separation  distance  s  and  another  on 
the  angle  y  between  the  centerline  0102  and  the  x  axis  (Fig. lb).  Let  A'jj  (i,j  =  1,2,3  4)  be  the 
added  masses  evaluated  in  the  polar  coordinates  corresponding  to  the  unit  velocities  U,,  where  U[ 
and  LA  are  velocities  of  body  1  and  body  2  respectively  along  the  direction  of  the  centerline  Oj02, 
U3  and  LV  are  the  corresponding  velocity  components  perpendicular  to  O1O2  (Fig. lb).  Since 
body  1  and  body  2  ire  symmetric  with  respect  to  the  centerline  O1O2,  a  sign  change  in  U'3  should 
not  change  the  fluid  kinetic  energy  in  terms  of  Ay.  Therefore  A't3=0.  Similarly,  A'14  =  A’23  = 
A  24  =0.  Based  on  the  invariance  of  the  kinetic  energy  to  the  coordinate  transformation,  we  obtair 
the  relation 

Ait  —  A'ucos2y  +  A'33  sin2y,  A33  =  A'nsin2y  +  A’33  cos 27, 

A;:  =  A  22Cos2y  +  A 44  sin2y,  A 12  =  A’i2Cos2y  +  A  34sin2y, 

A 13  =  (A'n  -  A'33)  sinycosy,  A23  =  (A' 12  -  A'34)siny  cosy.  (13) 


The  denvauves  of  Ay  with  respect  to  c  and  h  are  then  expressed  as  functions  of  s  and  y, 


dAp 

3c 


3Ai  1  0  3 A'33  .  _ 

[ - cos-y-t- - sin^-y]  cosy 

3s  3s 


2f-A'i  ]  +A'33)cosy  siny(-^~). 


3A[  1 

3h 


,3A' 


1 1 


3s 


i  3 A'33  .  ,  ,  .  „  . ,  . ,  cosy 

cos-y-t - sin-yl  siny  +  2(-A  1 1  +A  -<3)cosy  siny(— — )  ,  etc. 

3s  s 


n 


Relying  on  the  foregoing  transformation,  we  can  consider  the  added  masses  due  to  centroidal  and 
transversal  motions,  separately,  at  any  time  instant. 

III.  ADDED  MASSES  AND  THEIR  DERIVATIVES 

The  added  masses  of  two  bodies  moving  through  an  inviscid  and  incompressible  fluid  can 
be  determined  in  terms  of  hydrodynamic  singularities  within  each  body,  by  means  of  a 
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generalization  of  Taylor's  added-mass  theorem!9!.  For  a  pair  of  three-dimensional  solids  moving 
in  any  manner  except  pure  rotations,  the  generalized  Taylor’s  added-mass  formula  is  given  as 

Aiajp  +  5ap5ij  M  p  =  4rcp  [  J  -xj  XiapdV  +  I  (miap  x°j  +  (iiajp)  ] 

VP 

(no  sum  on  P),  (15) 

where  A;ajp  are  added  masses!7!  of  body  (3  along  the  j-th  direction  due  to  the  unit  motion  of  body 
a  in  the  i-th  direction,  M  p  is  the  mass  of  the  fluid  displaced  by  body  P,  A^ap  and  mjap  are  the 
volume -distributed  source  density  and  the  isolated  source  strength,  respectively,  inside  body  P  due 
to  the  i-th  velocity  component  of  body  a,  (ijajp  is  the  strength  of  an  isolated  dipole  in  the  j-th 
direction  inside  body  p  associated  with  the  i-th  velocity  component  of  body  a,  xj  and  x°j  are  the 
j-th  local  coordinates  of  A^p  and  rriiap  respectively  with  respect  to  body  P,  the  integration  is  over 
the  volume  of  body  P,  Vp,  and  the  summation  is  over  all  isolated  singularities  inside  Vp.  The 
factor  4  in  equation  (15)  should  be  changed  to  2  for  a  pair  of  two-dimensional  bodies.  The 
distributed  doublets  do  not  appear  explicitly  in  (15)  because,  by  integration  by  parts,  they  are 
equivalent  to  a  distribution  of  sources  in  Vp  plus  a  surface  distribution  of  sources  on  the 
boundary!10!.  For  the  present  two-dimensional  problem,  (15)  may  be  simplified  as 


(no  sum  on  j),  (16) 

where  the  corresponding  notations  are  defined  in  a  similar  manner  to  the  foregoing  ones. 

When  two  circular  cylinders  move  along  the  centerline  o\02,  with  velocities  U j  and  Ui 
respectively,  the  method  of  successive  images  (Hicks!5!,  Miillert8!)  can  be  applied  to  obtain  added 
masses.  This  method,  based  on  Taylor's  added-mass  formula  (16)  and  the  circle  theorem, 
produces  simple  expressions  for  added  masses  in  terms  of  the  strengths  of  isolated  doublets  inside 
circles  which  represent  an  equivalent  velocity  potential  due  to  the  motion  of  solids. 
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Consider  first  the  system  of  hydrodynamic  singularities  corresponding  to  the  unit  motion  of 
circle  1  along  the  centerline  0102,  U'i  =  1,  U'2  =  U3  =  U4  =  0  (Fig.  lb).  If  circle  2  were  absent, 
the  velocity  potential  could  be  represented  by  a  doublet  of  strength  |i0  =  a2  located  at  oj  in  the 
direction  of  U  j  (Fig. 2).  However,  the  presence  of  circle  2  violates  the  boundary  condition  on 
surface  2  and  needs  an  image  doublet  of  strength  ji!  =  -  a 2b2/s2  at  the  inverse  point  Pi  inside  circle 
2  to  satisfy  the  boundary  condition  on  circle  2,  where  the  minus  sign  for  |ii  indicates  that  the 
doublet  is  in  the  opposite  direction  of  U'i  and  the  point  P]  is  on  the  centerline  0102  with  P1O2  = 
b2/s  (Fig.2).  The  presence  of  this  image  doublet  violates  the  boundary  condition  on  surface  1  and 
requires  another  isolated  image  doublet  inside  circle  1  and  so  on.  The  general  expression  for  the 
strengths  and  positions  of  the  2n-th  (in  circle  1)  and  the  (2n+l)-th  (in  circle  2)  doublets  is  given  by 
the  iterative  formula  (Rouse  et  alJ11! ) 

Po  =  ^-o  =  S, 


P2n+t  —  '  P2n 


f  b2  h 


^•2n+l  =  S  - 


b2 

^•2n 


P2n+2 - P2n+1 


r  a2  ' 

^2n+l2y 


W2  =  s--J— ,  (n=0,l,2,...),  (17) 

^-2n+l 


where  A.2n  denotes  the  distance  from  02  to  the  n-th  inverse  point  in  circle  1  and  A.2n+i  that  from  oj 
to  the  (n+l)-th  inverse  point  in  circle  2,  all  inverse  points  lie  on  the  centerline  O1O2.  As  n  tends  to 
infinity,  the  system  of  doublets  generates  the  same  velocity  potential  exterior  to  the  circles  as  that 
due  to  the  unit  motion  of  circle  1  along  O1O2. 
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To  find  the  closed-form  expression  of  added  masses  from  the  iterative  formula  (17),  it  is 
necessary  to  represent  the  doublet  strength  at  the  n-th  iteration  in  terms  of  radii  a,  b,  and  the 
separation  distance  s.  A  series  expansion  of  added  masses  of  two  spheres  due  to  their  centroidal 
motion  was  derived  by  Hicksl5!  and  Herman!12!  by  means  of  the  method  of  continued  fraction. 
Their  result  was  also  generalized,  with  slight  modifications,  to  the  case  of  two  or  multiple  cylinders 
by  Miilleri8!  and  Yamamoto!3!.  However,  an  alternative  expansion  is  used  in  the  present  study  so 
that  we  may  find  the  derivatives  of  added  masses  with  respect  to  s  easily. 

Equation  (17)  can  be  reduced  to 


!A>n  =  a2 


(a  b)2n 

(X,0XlA.2--.X2n-l)2 


(18a) 


ti2n+l  — 


(a  b)2n+2 

(X0X1X2...X2n)2 


(n  =  0,  1,  2,  ...  )  . 


(18b) 


For  every  n,  the  location  of  the  inverse  point  satisfies  the  relations 
X2n^-2  n+1  -  sX2n  -  b2,  X.2  n+lX2n+2  -  a2  . 

After  some  manipulation,  it  can  be  found  that 

,  ,  ,  a 2 ( r i n  -  r2n)  +  (rtn+1  -  r2n+l) 

^-2n-1^2n-2  —  ~  ,  _  > 

ri  -  r2 


(19) 


(20a) 


and 


X2n^-2n-l — Xj  — 


rjn+1  _  f2n+l 

ri  -  r2  ’ 


(20b) 


where 

ri  2  =  ~  [  (s2  -  a2  -  b2)  ±  V  (s2  -  a2  -  b2)2  -  4  a2b2  ]  . 
2 


(21) 
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The  n-th  image-doublet  strength  inside  circle  1  is  then 


_ (a  b)2n _ 

[a2(rin'1+rin'2r2+...+r2n'1)+(rin+rin'1r2+...+r2n)]2 


(22a) 


and  that  inside  circle  2  is 

_  _ (a  b)2n+2 _ 

^2n+1  s2(rin+rin-lr2+ ---+r2n)2 


(22b) 


For  a  fixed  n,  the  limits  of  sequences  {q.2n) .  {|i2n+i }  as  s  approaches  (a+b)  are  derived  as 


a2  b2 

lim  (i2n  = - r  ,  (23a) 

s—»(a+b)  (n(a  +  b)  +  b) 


a2  b2 

lim  42n+l  =  - - V5 - 2  •  (23b) 

s-»(a+b)  (a+b;(n+l)z 


(p.2n)  and  ((i;n+i }  are  uniformly  convergent  for  any  values  of  s  in  the  region  [a+b,°°)  since  in  the 
limit,  s — >(a+o),  (23a,b)  decay  by  the  rate  of  1/n2  and  for  s  >  (a+b), 

^  =  f<l  and  f<l. 

M-2n  V^2n^2n+lJ  M-2n-l  V^2n-1^2nJ 


Substitution  of  P2n  and  p.2n+i  into  Taylor's  added-mass  formula  (16)  yields  the  expression 

of  added  masses  due  to  the  centroidal  motion  of  two  circles  in  the  relative  polar  coordinates 

00  2 

A'n  =  2np[  Sp2n-  y]  , 
n=0 


A'i2  =  2rcp  X  M-2n+l  , 
n=0 
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A'22  =  2  Tip  [ 


00  b2 
IH*2n-yl, 

n=0 


A '21  =  A’i2 


(24) 


where  p  is  the  density  of  the  fluid,  |i*2n  denotes  the  doublet  strength  inside  circle  2  due  to  the  unit 
motion  of  itself  along  the  centerline  0102,  U'2=l,  U’i=U'3=U4=0,  and  is  obtained  directly  by 
interchanging  a  and  b  in  equation  (17).  The  limiting  values  of  added  masses,  as  s  approaches  (a  + 
b),  are  obtained  from  expressions  (23a, b)  and  (24) 


r  a2b2  “ 
hm  A  n  =  2  tt  p  [ - r£ 

s— »(a+b)  (a  +  b)  n=0 


(25a) 


2  U2  2  k2 

lim  A'12  =  -  2  tt  p  C  (2)  =  -  3.28968  Tt  p  J— ^  ,  (25b) 

s— »(a+b)  (a+b)  (a+b) 


where  £(2)  is  the  zeta  function.  From  (22a,b)  and  (24),  we  note  that,  as  s  increases  from  (a+b)  to 
infinity,  the  added  mass  A'n  decreases  monotonically  from  its  limiting  value  given  in  (25a)  to  a 
constant  7tpa2  and  the  interaction  added  mass  A'12  increases  monotonically  from  the  value  given  in 
(25b)  to  zero. 

To  evaluate  the  derivatives  of  added  masses  directly  from  equation  (24)  by  differentiating 
each  term  of  sequences  {(i2n)  anti  {(i2n+i)  with  respect  to  s,  we  need  to  prove  that  in  a  certain 
region  of  s,  the  sequences  {dji2n/ds}  and  {dp.2n+i/ds)  are  uniformly  convergent.  By  taking  the 
derivatives  of  p.2n  and  P2n+i  directly  from  (18),  we  have 


tiM-2n 

ds 


a2  a2n  b2n 


(^2n-1^2n-2----^o)3 


(^*2n-l ^-2n-2 — ^o)* 


(26a) 
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and 


dM-2n+l 

ds 


2  (a  b)2n+2 


[1+- 


S3(^2n^2n-l--”A.l)2  ^2n^2n-l----^l 


^  (X2n^2n-l----^l)]  •  (26b) 


Let  n0  be  the  product  A.2n-i^2n-2---^o.  Since ,  from  (21), 


dri  __  2s  ri  dr2  _  2s  r2 
ds  ~  n-r2  ’  ds  ~  n-r2  ’ 

the  differential  part  of  (26a)  becomes 


and  rir2=(ab)2, 


(27) 


dllp  3n0  dri  ailp  dr2 
ds  ~  3n  ds  +  5r2  ds" 


2s 

(n-r2)3 


n-I  n-1  „  n+I  n+I  n+2  n+2  n  n 

{-(ab)2a2(n+l)(r  -r  )+a-(n-l)(r  -r  )+n(r  -r  )-(n+2)(ab)2(r -r  ) } . 

12  12  12  12 


The  limiting  value  of  (26a)  for  fixed  s  >  (a+b)  is 


(28) 


Um  dM-2n  =  _  iim  4  a2(ab)2n  s  n 
n_^co  ds  n->~  n02(n-r2) 

Based  on  the  ratio  test. 


=  0. 


(29) 


lim  dM-2n+2/ds _ (ab)2  < 

n~*°°  c4Wds  "a2nX2n+i)2  <  ’ 


we  conclude  that  the  sequence  {d(j.2n/ds}  is  uniformly  convergent  in  the  region  (a+b,  <»). 
Similarly,  we  can  prove  that  the  sequence  {dp.2n+i/ds}  is  also  uniformly  convergent  in  the  same 
region.  Thus,  the  derivatives  of  added  masses  (24)  with  respect  to  s  are  equal  to  the  summations 
of  derivatives  of  the  corresponding  doublet  strengths  for  (a+b)  <  s  <  <». 

From  (17)  and  (24),  we  have 

_  a  d^o  , 

ds  ~u’  ds  ~ 
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d42n+L=  d)i2n  f  b2 
dS  '  dS 


\2n2Jb2  Yd^ 


J 


x2n3  K  ds  r 


d^2n+l  _  x  +  (  b2  YdX2n 


ds 


V 


A-2n2 


ds 


dM-2n+2 

ds 


d}l2n+l 

f  a2 

ds 

^2n+l2 

dX2n+2 

ds 


/  a2  \ 

f dA^n+A 

^2n+l2y 

l  ds  J 

(30) 


and 


dA'n 

ds 


=  2irp 


”dp2n 

Ln  ds 
n=0 


dA' 


12 


ds 


=  2?tp  I 


dp2n+l 


n=0 


ds 


dA'22 

ds 


=  27tp 


dM-*2n+l 


n=0 


ds 


dA'2i_dA'12 
ds  ds 


(31) 


The  region  of  convergence  for  derivatives  of  added  masses  in  (31)  is  open  on  the  left  since  the 
limiting  value  of  (29)  is  not  zero  as  s  approaches  (a+b). 

The  above  analysis  for  added  masses  and  their  derivatives  indicates  that,  as  the  separation 
distance  s  tends  ro  (a+b),  the  kinetic  energy  of  the  fluid  due  to  the  motion  of  solids  is  finite  but  the 
hydrodynamic  interaction  forces,  which  are  functions  of  the  derivatives  of  added  masses  with 
respect  to  s,  approach  infinity. 

When  two  circles  make  transversal  motions  perpendicular  to  the  centerline  0i02,  the 
strengths  and  locations  of  doublets  in  each  circle  can  be  determined  in  the  same  way  as  that  used 
for  the  centroidal  motion.  Thus,  we  have 

A'n  =  A'33,  A'22  =  A'44,  A'i2  =  -A'34  .  (32) 
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Since  the  summations  of  doublet  strengths  and  their  derivatives  with  respect  to  s,  given  by 
iterative  formulas  (17)  and  (30),  are  convergent  to  the  corresponding  solutions  of  added  masses 
and  their  derivatives,  we  can  calculate  them  numerically  based  on  these  iterative  formulas  and 
obtain  the  results  as  accurate  as  we  want  by  increasing  the  index  n  to  a  sufficiently  large  number. 
Practically,  if  s  is  larger  than  1.5(a+b),  an  accuracy  of  five  significant  figures  can  be  achieved 
within  a  few  iterations  and  for  (a+b)  <  s  <  1.01  (a+b),  the  same  accuracy  may  need  several 
hundred  iterations. 

IV.  NUMERICAL  SOLUTION  OF  ADDED  MASSES  AND  THEIR  DERIVATIVES 
The  application  of  the  afore-mentioned  analysis  is  limited  by  the  body  geometry  and  flow 
conditions  since,  in  general,  the  velocity  potential  cannot  be  expressed  in  terms  of  isolated  interior 
singularities.  For  bodies  of  various  geometries,  however,  we  may  have  to  evaluate  added  masses 
and  their  derivatives  numerically  in  terms  of  certain  distributed  singularities  either  inside  bodies  or 
on  the  body  surfaces.  Landweber  &  Chwangf7!  has  developed  a  boundary-integral  method  for 
two-body  interaction  problems  based  on  the  generalized  Taylor's  formula  and  the  fundamental 
relation  between  the  velocity  potential  and  normal  velocities  on  solid  boundaries.  The  accuracy  of 
their  method  depends  on  the  simultaneous  solutions  of  a  set  of  Fredholm  integral  equations  of  the 
second  kind.  Due  to  the  singular  behavior  of  the  kernels  in  these  equations,  some  numerical 
inaccuracy  in  the  solutions  may  be  significant  when  two  bodies  are  very  close  to  each  other.  In 
order  to  investigate  the  reliability  and  to  improve  the  accuracy  of  numerical  results,  we  shall 
consider  a  pair  of  circles  again  and  compare  the  numerical  results  with  exact  solutions. 

From  the  work  of  Landweber  &  Chwangl7!,  added  masses  due  to  the  general  translational 
motion  of  a  pair  of  two-dimensional  bodies  are  given  by 

A,aj(3  =  27tp  j  XjajapdSp  -M'jjSijSap  (no  sum  on  (5),  (33) 
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where  p  is  the  density  of  the  fluid,  M'p  ((3=1,2)  is  the  mass  of  the  fluid  displaced  by  body  (3,  ajap 
denotes  the  surface  source  distribution  on  body  (3  due  to  the  i-th  unit  velocity  component  of  body 
a,  and  Xj  is  the  coordinate  of  the  surface  source  on  body  [3.  For  the  relative  motion  of  two  circles 
concerned  in  the  present  study,  the  unknown  surface  source  distributions  satisfy  a  set  of  integral 
equations 


it  E(P)  + 


r-/™  C0S  Vl2  TT'  dx  .  T  T-  dy 

F(Q)  — 5 - -dS2  =  Uir —  +  U3— — 

K2i  dm  dni 


(34a) 


JtF(Q) 


cos  Wi 1  .  dx  ,  dy 

E(P')  -—^dSi  =U2— +U4-^ 
k12  dn2  dn2 


(34b) 


where  E(P)  denotes  the  unknown  source  strength  at  point  P  on  circle  1,  F(Q)  denotes  that  at  point 
Q  on  circle  2,  is  the  distance  between  the  source  point  on  body  i  and  the  field  point  on  the 
surface  of  body  j,  n,  indicates  the  outward  normal  direction  at  the  field  point  on  the  surface  of  body 
i,  and  \jqj  is  the  angle  between  nj  and  Rji  (Fig. 3a).  In  the  derivation  of  equation(34a,b),  integrals 
involving  the  kernel  (cos  yn)/Rn  or  (cos  \|/22)/R22  reduce  to  zero  for  circles,  since  they  are 
constant  on  circles  1  and  2  and  the  integrals  of  E  and  F  over  closed  bodies  vanish  identically. 

To  obtain  the  added  masses,  we  first  set  U  1  =  1,  U’2  =  U  3  =  U'4  =  0.  The  solution  of 
(34)  for  this  condition  yields  the  source  distributions  for  evaluating  A’n  and  A\2, 

2k 

A'i  1  =  2^pa2  J  E(a)  cosa  da  -  prca2  ,  (35a) 

o 

2n 

A  12  =  2npb2  j  F(9)  cos9  d9  .  (35b) 

o 


Similarly,  A  22  is  obtained  by  using  the  source  density  F  determined  from  the  condition  that  U'2  = 
1  and  U'i  =  U3  =  U4  =  0, 
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2k 

A 22  =  2irpb2  J  F(0)  cos0  d0  -  pttb2.  (35c) 

o 


From  (32)  and  the  symmetry  argument  leading  to  (13),  we  have 


> 

II 

> 

T 

II 

A'23  =  A  24  =  0, 

(35d) 

A  ii  =  A33, 

A'22  =  A44,  A  12  =  -  A  34. 

(35e) 

The  dimensionless  added-mass  coefficients  may  be  defined  as 

ky  =  A'ij/  M'i.  (36) 

All  integrations  in  equations  (34a, b)  are  regular  for  s  >  (a+b).  However,  by  comparing 
the  analytical  solution  of  added  masses  from  equations  (17)  and  (24)  with  the  numerical  result  from 
integral  equations  (34)  and  (35)  (Table  1,  columns  with  superscripts  "e"  and  "a"),  we  recognize 
that  a  considerable  deviation  in  the  numerical  result  takes  place  when  s  =  a+b+5,  5«  1.  The 
numerical  error  comes  mainly  from  two  sources,  one  due  to  the  ill-conditioned  kernels  which  have 
a  large  and  steep  peak  value  of  order  o(l/5)  in  the  region  around  a=iz  and  0=0  (Fig. 3a),  and 
another  due  to  the  discretization  of  the  integral  equations  into  a  set  of  linear  algebraic  equations.  In 
order  to  remove  the  peaks,  we  subtract  a  function,  which  has  the  same  peak  value  and  similar 
behavior  around  a=7t  and  0=0,  from  the  integral  and  then  add  the  exact  integration  of  the  function 
back  to  the  equation.  Thus,  equation(34a,b)  is  modified  to 

2k  2k 

nE(a)+b  f  [F(0’)  C°^-¥l-2  (a,0')  -  F(0)  gl(a,0’)]d0’+b  J  F(0)  gl(a,0')d0' 

O  K21  O 

=  U'icosa+U'3sina,  (37a) 

2K 

7tF(0)+a  j  (E(a’)  -°p^21  (a',0)  -  E(tc)  g2(a’,0)]da’  +  a  J  E(n)  g2(a’,0)d<x' 
o  Kl2  o 
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=  U^cosO+U^sinG 


(37b) 


where 


cos  V/12 
R21 


a  -  bcos(0'  -  a)  +  s  cosa 


2  2  ' 
(s  -  b  cos0'  +  a  cosa)  +  (b  sin0'  -  a  sina  ) 


(38a) 


COS  V21 

Rl2 


b  -  acos(0  -  a')  -  s  cos0 


2  2 
(s  -  b  cos0  +  a  cosa')  +  (b  sin0  -  a  sina'  ) 


(38b) 


For  simple  kernels  like  (38a,b),  we  may  let  gi  and  g2  be  (cos  ¥i2)/R2i  and  (cos  V2O/R12 
respectively.  Then  integrations  of  gi  and  g2  on  circles  2  and  1,  respectively,  give  (Fig. 3) 


271 


where 


and 


b  /  gi(a,0')d0’=  ~  w-b  cos(a  -  p), 
o  u 

2?r  2tt  a 

a j  g2(a',0)da'= -—r^cos(0  +  p*), 
o  u 


d2  =  s2  +  a2  +  2  a  s  cosa,  d*2  =  s2  +  b2  -2  s  b  cos0, 


P=±cosVii^_aiy  P.=±co,>M^ 


(39a) 

(39b) 

(40) 

(41) 


The  positive  branches  of  p  and  p*  correspond  to  0  <  a  <  7t  and  0  <  0  <  k  respectively,  and  the 
negative  ones  to  n  <  a  <  2k  and  7t  <  0  <  2k.  Numerical  results  of  ku  and  k^,  defined  by 
equation  (36),  based  on  the  modification  are  listed  in  Table  1.  From  columns  with  superscript  "c" 
in  Table  1,  we  note  immediately  the  significant  improvement  based  on  equation  (37). 

Although  the  integrands  in  equation  (37a,b)  vanish  at  0=0  and  a=7t,  they  have  not  been 
quite  smoothed  in  that  neighborhood.  In  this  peak  region,  the  difference  of  source  strengths 
F(0’)-F(O)  or  E(a’)-E(7t)  is  not  large,  but  its  product  with  the  ill-conditioned  kernel  can  become 
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very  large,  especially  for  the  centroidal  motion,  such  that  it  affects  the  evaluation  of  the  integrals 
significantly.  The  accuracy  may  be  further  improved  by  choosing  an  appropriate  quadrature 
formula.  It  has  been  proved  by  Atkinson^13'  that  if  the  integrand  is  a  periodic  function  and  the 
domain  of  integration  is  evenly  discretized,  the  trapezoidal  rule  provides  the  best  numerical 
estimation  of  the  integral.  However,  it  is  clear  that  we  should  put  more  points  in  the  peak  region 
and  less  points  in  the  region  where  the  integrand  slopes  gently.  In  the  present  calculation,  we 
subdivide  the  region  of  integration  into  two  parts,  one  contains  the  sharp  peak  and  the  other  covers 
the  rest(Fig.4),  and  in  each  subdivision,  evaluate  integrals  in  Eq.  (37a, b)  by  the  Gaussian 
quadrature  formula,  in  order  to  compare  numerical  results  of  the  integration  in  the  peak  region, 
obtained  by  the  Gaussian  quadrature  formula  and  by  the  trapezoidal  rule,  with  the  exact  solution, 
we  examine  the  integral 


in 


I  = 


c^Vi2,«,e  >  de , 

^21 


(42) 


where  (cos  Vi2)/R-2i  is  given  by  (38a)  and  a  is  the  angle  between  01O2  and  0]P(Fig.3b).  At  s  = 
a+b+5  and  a/b  =  1.0,  the  integrand  in  (42)  has  a  peak  of  the  order  of  o(l/5)  around  a  =  re  and  9  = 
0.  The  results  are  given  in  Fig.  5.  From  this  figure,  we  observe  that  the  deviation  of  the 
numerical  results  obtained  by  the  Gaussian  quadrature  formula  from  the  exact  solution  is  much 
smaller  than  that  by  the  trapezoidal  rule.  For  a/b  =  1  and  various  values  of  s,  the  numerical  results 
of  added-mass  coefficients,  kn  and  ki2,  obtained  by  the  trapezoidal  rule  with  80  points  and  by  the 
Gaussian  quadrature  formula  with  60  points  along  each  circle,  are  listed  in  Table  1  (columns  "b" 
and  "c”).  The  latter  has  only  less  than  1%  error  at  s=2.01,  while  the  former  does  not  converge  at  s 
=  2.01  and  has  1.4%  error  for  kj2  at  s=2.04.  The  size  of  the  subdivision  on  each  circle  depends 
surely  on  the  separation  distance  s,  since  the  peak  becomes  sharper  as  s  is  closer  to  (a+b).  In  the 
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present  calculation,  the  divisions  are  roughly  estimated  and  their  sizes,  measured  by  the  angle  0) 
for  a=b,  are  listed  in  Table  1. 

In  order  to  find  derivatives  of  added  masses  with  respect  to  the  separation  distance  s  in 
terms  of  the  boundary  integrals,  it  is  necessary  to  evaluate  the  derivative  of  the  surface  source 
distribution  with  respect  to  s.  For  two  circles,  the  derivatives  are  given  implicitly  by  a  pair  of 
Fredholm  integral  equations  of  the  second  kind, 

rcEs(a)+b  J  [ Fs(0  )  -  FS(Q) ] C°p-Xj/ 1  -  (a.0’)  d0’  + 
o  ~i 


2k 


b  j  [ F(0 ')  -  F(0)]  ^ 

o  a\ 


COS  VJ/ 1  2 


R 
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(a,0')  a0'  =  -  P(a), 


COS  tit'i  I 

TtFsi.0)+aj  [Es(a')  -  Esctt) J — - — —  (a',0)  da'  + 
o 


(43a) 


In 


a  j  |E(a  )  -  E(rt)| 
o 


d  / cos  V:  j 
<*\  Ri: 


(a',0) 


da'  =  -  Q(9) 


(43b) 


where  Eqai  and  FqO)  stand  for  dE(a)/ds  and  dF(0)/ds  on  circles  1  and  2  respectively,  and 


In 


P(a)  =  b  {  [  FdO)C°p  Xl/—  (a,0')  +  F(0)jj ^ 12  (a,0’j  |  1  d0' 

q  Kij  uM  ^21 


—j  { Fs(0)  cos(a  -  (3)  -  [a  cos(3  +  s  cos(a-P)] } . 
a  d2 


(44a) 


COS  \p2  1 


df COS  \|/21 


Q(0)  =  a  j  i  Es(n)— 5  1  -  (a',9)  +  EfnVm  ■  D-  -  (a',0)  ]  da' 

n  k12  asl  Ki2  I 


2  k  a 


E(jt) 


jlj —  {  Es(7C)  cos(0  +  p*)  +  -j— [b  COsP*  -  s  cos(0  +  P 


■)]}. 


(44b) 
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Following  the  procedure  in  solving  (37a,b),  we  subdivide  the  region  of  integrations  into  two  parts, 
and  in  each  subdivision,  evaluate  integrals  in  (44a,b)  by  the  Gaussian  quadrature  formula.  The  set 
of  integral  equations  in  (44a, b)  is  then  solved  by  the  Gauss-Seidel  iterative  method  After  the 
determination  of  functions  Es(a)  and  Fs(0),  the  derivatives  of  added  masses  can  be  obtained  by 
simple  integrations.  For  example,  from  (35a)  and  (35b),  we  have 

dA'  2n 

— t  1 1  =  2rcpa2  f  Es(ot)  cosa  da  , 
as  o 

dA'  lK 

■~t  12  =  2ttpb2  J  Fs(0)  cos0  d0  .  (45) 

as  o 

An  accurate  solution  of  (43a, b)  is  more  difficult  to  obtain  than  that  of  (37a, b)  when  two 
bodies  are  close  to  each  other,  since  the  kernels  d(cos\j/i2/R2i)/ds  and  d(cos\j/2i/Ri2)/ds  in  (44a, b) 
have  peaks  of  the  order  of  o(l/82)  in  the  region  around  a=K  and  0=0.  To  solve  this  problem,  we 
simply  add  more  points  in  the  peak  region  when  s  is  close  to  (a+b).  For  a/b=1.0,  all  integrals  in 
equations  (44a, b)  are  evaluated  by  the  160-point  Gaussian  quadrature  for  2.0  <  s  <  2.1  and  by  the 
60-point  Gaussian  quadrature  for  s  >  2.1.  The  numerical  values  of  the  derivatives  of  added-mass 
coefficients  are  compared  with  the  exact  ones  in  Table  2.  At  s  =  2.01,  the  error  for  the  derivative 
of  kj  [  with  respect  to  s  is  2.6%  and  that  for  dki^/ds  is  2.5%.  For  s  >  2.02,  however,  the  results 
>f  (45)  based  on  solutions  of  (43a, b)  agree  with  the  exact  values  up  to  j  significant  figures. 
Therefore,  they  may  be  regarded  as  exact  for  all  practical  purposes. 

V.  NUMERICAL  RESULTS  ON  TRAJECTORIES 

For  a  pair  of  circular  cylinders,  the  added  masses  and  their  derivatives  with  respect  to  s 
have  been  evaluated  exactly  by  the  method  of  successive  images  and  approximately  by  the 
boundary  integral  method.  By  (13)  and  (14),  the  results  obtained  in  the  relative  polar  coordinates 
can  be  transformed  into  the  Cartesian  coordinates  in  which  the  motions  of  bodies  are  described. 
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With  A13  =  0  and  An  =  A33,  the  Lagrange’s  equations  of  motion(l  la,b)  are  decoupled  and 
simplified.  The  instantaneous  velocity  components  u  and  v  of  circle  1,  which  are  functions  of 
added  masses  and  their  derivatives  with  respect  to  c  and  h(Fig.  1),  are  given  by  the  first  integration 
of  (1  la,b)  with  respect  to  t.  The  trajectories  of  the  moving  circle  around  the  fixed  one,  are  then 
integrated  from  u  and  v.  In  the  calculation,  the  origin  of  the  absolute  coordinate  system  is  located 
at  the  center  of  the  fixed  body.  We  assume  that  the  uniform  flow  conveying  the  drifting  body  is 
U0=  1.0  along  the  x  direction  at  infinity.  When  two  circles  are  far  away  from  each  other,  we  may 
specify  the  initial  condition  used  in  the  solution  as  u0  =  1.0,  v0  =  0.0  for  s0/(a+b)  >  10. 
Numerically,  we  discretize  the  time  t  into  a  number  of  intervals,  St;  (i=0, 1,2, ... ),  and  assume  that 
within  each  time  interval,  velocities  u  and  v  are  constants  which  are  the  same  as  those  at  the 
beginning  of  the  interval.  The  lengths  of  St  are  taken  as  St;  =  0.1  when  x  <  -3  and  St;  =  0.01 
when  x  is  larger  than  -3. 

The  computation  is  started  at  the  initial  position  (x0,  y0),  where  xG  =  -20  for  b  =  1,  a  <  1 
and  y0  varies.  At  the  beginning  of  each  interval  St;,  u;  and  v;  are  obtained  by  the  fourth-  order 
Runge-Kutta  integration  of  (1  la,b)  in  terms  of  velocities  in  the  previous  interval,  u;_i  and  Vj.j,  the 
added  masses,  and  their  derivatives.  By  the  end  of  this  interval,  the  new  position  of  the  drifting 
circle  is  given  by 

X!  =  X;_1  +  UjSt;  and  y;  =  y;.;  +  VjStj  . 

This  process  is  repeated  until  the  separation  distance  s  is  equal  to  or  less  than  (a  +  b),  at  which 
impact  between  two  bodies  takes  place.  The  radii  a  and  b  are  normalized  by  that  of  the  fixed 
cylinder.  Various  ratios  of  a/b,  depending  on  different  physical  problems,  are  considered  in  the 
calculation. 

We  first  consider  a  circular  cylinder  of  ice  with  radius  a  =  0. 1  and  density  p=0.91 ,  moving 
around  a  fixed  cylinder  of  radius  b=1.0.  At  yo=0. 1,  0.3,  etc.,  the  trajectories  of  the  ice  particle 
are  plotted  in  Figures  6a  to  6e  for  different  density  ratios  of  the  body  to  the  fluid  medium.  Fig. 6a 
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corresponds  to  the  case  of  an  ice  particle  moving  in  fresh  water  and  Fig.6e  shows  the  same  ice 
particle  in  an  air  flow.  The  ratios  2.0,  5.0  and  10.0  correspond  to  various  fluid  media  and 
illustrate  the  change  of  trajectories  with  respect  to  the  density  ratio.  We  note  from  these  figures 
that  the  lower  the  density  ratio  is,  the  larger  the  curvature  of  trajectories  would  be.  In  the  case  of 
an  ice  particle  carried  by  an  air  flow,  its  trajectories  are  almost  straight  lines  which  differ  from  the 
streamlines  greatly.  It  means  that  in  a  fluid  of  high  density,  the  hydrodynamic  interaction  due  to 
the  presence  of  a  second  body  has  a  large  effect  on  the  body’s  motion. 

In  some  physical  problems,  for  example,  in  the  ice  coating  process,  it  is  important  to 
determine  whether  or  not  a  drifting  body,  conveyed  by  a  current,  can  impact  with  a  fixed  body. 
This  physical  property  can  be  expressed  by  a  "collection  coefficient"  E  which  is  defined  as  the  ratio 
of  a  critical  initial  position  y*0,  below  which  the  drifting  body  will  impact  with  the  fixed  body,  to 
the  radius  of  the  fixed  body.  y*0  is  surely  depending  on  the  ratios  of  a/b  and  pa/p.  In  this  type  of 
problems,  the  size  of  the  moving  body  is  much  smaller  than  that  of  the  fixed  one.  Fig.7  shows 
the  result  of  the  collection  coefficient  for  a  pair  of  cylinders,  0.01  <  a/b  <  1,  in  fluids  of  different 
densities.  Moving  bodies  of  very  small  size  are  not  considered  since  the  Reynold's  number 
becomes  quite  small  such  that  the  inviscid-fluid  theory  is  not  applicable. 

Regarding  to  a  floating  body  in  sea  water,  we  can  consider  the  interaction  problem  between 
an  ice  floe,  idealized  as  a  cylinder,  and  a  fixed  cylindrical  offshore-structure.  Fig. 8  shows  the 
trajectories  of  a  circular  cylinder  in  sea  water.  Here  the  density  ratio  pa/p  is  fixed  at  0.89,  which 
corresponds  to  a  floating  ice  in  sea  water.  The  radius  ratio  a/b  is  0.5  in  Fig. 8a  and  1.0  in  Fig. 8b. 
We  note  that  the  trajectory  of  a  large  body  is  quite  flat  in  comparison  with  the  one  of  a  small  body. 
It  is  understood  that  the  inertia  effect  of  a  large  body,  which  prevents  the  body  getting  off  the 
straight  path,  predominates  over  the  hydrodynamic  interaction  force  due  to  the  presence  of  a 
second  body. 

For  a  variety  of  radius  ratios  a/b  from  0. 1  to  1 .0  and  a  fixed  density  ratio  0.89,  the  velocity 
components  u  and  v  of  the  moving  cylinder  are  shown  in  Fig.9.  The  initial  condition  used  in  the 
calculation  is  xQ  =  -20,  y0  =  0.5,  uo=1.0,  and  v0  =  0.0.  This  plot  is  consistent  with  the  physical 
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interpretation.  The  velocity  a  is  reduced  as  the  moving  body  approaches  the  fixed  body,  and  after 
a  certain  position,  it  increases  and  reaches  a  maximum  value  around  the  top  of  the  stationary 
cylinder.  On  the  other  hand,  the  velocity  v  increases  from  zero  to  a  maximum  value  in 
approaching  and  decreases  to  zero  after  the  cylinder  passes  over  the  fixed  one.  We  also  note  from 
Fig. 9  that  for  x  <  -2.5,  velocity  components  for  cylinders  of  different  radii  are  of  the  same  value 
practically.  It  indicates  that  for  two  bodies  apart  from  each  other  by  a  large  distance,  the 
hydrodynamic  interaction  does  not  have  any  effects  on  their  motion. 

VI.  CONCLUSIONS 

The  hydrodynamic  interaction  between  two  circular  cylinders  of  different  radii  moving 
relative  to  each  other  in  a  general  manner  is  investigated  based  on  the  generalized  Taylor's  added- 
mass  formula  in  terms  of  isolated  interior  doublets  or  surface  source  distributions.  The  velocity 
potential  of  the  fluid  is  expressed  in  terms  of  three  independent  added  masses  and  their  derivatives 
with  respect  to  the  separation  distance,  which  are  functions  of  relative  sizes  of  the  two  cylinders 
and  the  separation  distance.  By  extending  Hicks'  analysis,  series  expressions  have  been 
developed  for  derivatives  of  added  masses  based  on  the  successive  image  method.  In  the  general 
case  of  an  oblique  impact,  it  has  been  found  that  the  numerical  evaluation  of  the  derivatives  of 
added  masses  is  more  difficult  than  that  of  added  masses  themselves,  because  of  the  sensitive 
dependence  of  integrals  in  the  integral  equations  on  the  separation  distance. 

The  Lagrange's  equations  of  motion  governing  the  motion  of  a  drifting  cylinder  relative  to  a 
stationary  one  are  integrated  to  obtain  the  trajectory  of  the  moving  cylinder.  The  numerical  results 
exhibit  the  dependence  of  the  hydrodynamic  interaction  force  on  the  direction  of  the  flow  with 
respect  to  the  centerline  joining  the  centers  of  two  cylinders.  The  velocity  component  along  the 
centerline  produces  a  repulsive  interaction  force,  which  prevents  the  collision  of  two  bodies,  while 
the  component  perpendicular  to  it  produces  an  attractive  interaction  force.  Thus,  the  collision  or 
non-collision  prediction  for  general  oblique  motions  of  two  cylinders  is  much  more  complicated 
than  that  for  the  central-impact  problem. 
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Due  to  the  limitation  of  the  potential-flow  analysis,  the  afore-mentioned  results  are  only 
applicable  to  the  situation  where  the  effects  of  fluid  inertia  and  nonuniformity  of  the  flow  due  to  the 
presence  of  a  second  body  are  dominant. 
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Table  1.  The  Analytic  and  Numerical  Results  of  Added-Mass 
Coefficients  for  Two  Circles  fa/b  =  1 .0  ) 


s 

kne 

-k12e 

kna 

-ki2a 

kub 

-ki2b 

k  i  ic 

i 

7C 

ro 

o 

2.01 

1.3753 

0.7283 

*  *  *  * 

*  *  *  * 

**  ** 

**  ** 

1.3713 

0.7214 

12.0 

2.02 

1.3416 

0.6924 

**  ** 

**  ** 

1.1650 

0.5159 

1.3400 

0.6901 

12.5 

2.03 

1.3174 

0.6661 

1.8778 

1.2302 

1.2843 

0.6330 

1.3169 

0.6653 

13.0 

2.04 

1.2980 

0.6446 

1.3963 

0.7446 

1.2893 

0.6358 

1.2979 

0.6443 

13.5 

2.05 

1.2818 

0.6262 

1.3148 

0.6600 

1.2791 

0.6236 

1.2817 

0.6261 

14.0 

2.06 

1.2677 

0.6100 

1.2805 

0.6232 

1.2668 

0.6092 

1.2676 

0.6100 

14.5 

2.07 

1.2552 

0.5955 

1.2604 

0.6009 

1.2549 

0.5952 

1.2552 

0.5955 

15.0 

2.08 

1.2440 

0.5822 

1.2462 

0.5845 

1.2439 

0.5821 

1.2440 

0.5822 

15.0 

2.09 

1.2339 

0.5700 

1.2384 

0.5710 

1.2338 

0.5700 

1.2339 

0.5700 

15.5 

2.10 

1.2246 

0.5587 

1.2250 

0.5592 

1.2246 

0.5587 

1.2246 

0.5587 

16.0 

10.0 

m 

ImI 

KfflHtU] 

Bn 

iBI 

90. Q 

superscript:  e  -  exact  solution  from  the  method  of  successive  images 

a  -  numerical  solution  of  Eq(34)  by  the  trapezoidal  rule  with  80  points 
b  -  numerical  solution  of  Eq(37)  by  the  trapezoidal  rule  with  80  points 
c  -  numerical  solution  of  Eq(37)  by  the  Gaussian  quadrature  with  60  points 
d  -  (0  is  measured  by  degrees 

**.**  -  solution  not  convergent  within  50  iterations 
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Tabic  2.  Tb'31  Analvfc  and  Numerical  Revolts  of  the  Deri-. 

’stives 

of  Added-mass  Coefficients  for  Two  Circles  (a/b=1.0  and  w=20°) 

dk„ 

dkn 

dki2| 

dki2 

s 

ds  e 

ds  n 

ds  le 

ds  n 

2.01 

-4.22048 

-4.10983 

4.43434 

4.32371 

2.02 

-2.77334 

-2.77325 

2.98604 

2.98597 

2.03 

-2.13769 

-2.13769 

2.34932 

2.34931 

2.04 

-1.76190 

-1.76190 

1.97240 

1.97240 

2.05 

-1.50575 

-1.50575 

1.71691 

1.71691 

2.06 

-1.32129 

-1.32129 

1.52955 

1.52955 

2.07 

-1.17771 

-1.17771 

1.38485 

1.38485 

2.08 

-1.06291 

-1.06291 

1.26892 

1.26892 

2.09 

-0.96858 

-0.96858 

1.17347 

1.17347 

2.10 

-0.88942 

-0.88942 

1.09319 

1.09319 

Subscripts:  e  -  exact  solution  from  the  iterative  formulas  (30)  and  (31) 

n  -  numerical  solution  of  (43a,b)  and  (45)  by  the  Gaussian 
quadrature  formula  with  160  points 
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Figure  4.  Definition  of  regions  containing  peaks. 
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Figure  5.  Comparison  of  numerical  results  of  integral  I  in  equation  (42)  obtained 
by  the  trapezoidal  rule  with  80  points  and  by  the  Gaussian  quadrature 
with  60  points. 
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Figure  9.  Velocity  components  u  and  v  of  a  moving  cylinder  with  various  radius 
ratios. 
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